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Abstract 

We have observed bottom-charm mesons B c via the decay mode B^r — ► 
Jj^f^v in 1.8 TeV pp collisions using the CDF detector at the Fermilab 
Tevatron. A fit of background and signal contributions to the J/ip£ mass 
distribution yielded 2O.4+55 events from B c mesons. A fit to the same dis- 
tribution with background alone was rejected at the level of 4.8 standard 
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deviations. We measured the Bj mass to be 6.40 ± 0.39 ± 0.13 GeV/c 2 and 
the lifetime to be 0.46^q \q (stat.) ± 0.03 (syst.) ps. We measured the pro- 
duction cross section times branching ratio for B+ — > J ftp £ + v relative to that 
for B+ -► J/4>K+ to be 0.132±g;gf^ (stat.) ± 0.031 (syst.) tof 2 l (lifetime). 
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The .6+ meson is the lowest-mass bound state of a family of quarkonium states containing 
a charm quark and a bottom anti-quark [l]. Since this pseudoscalar ground state has non- 
zero flavor, it has no strong or electromagnetic decay channels. It is the last such meson 
predicted by the Standard Model. Its weak decay is expected to yield a large branching 
fraction to final states containing a. J ftp Non-relativistic potential models have been 

used to predict a B c mass in the range 6.2-6.3 GeV/c 2 |||7]|. In these models, the c and b 
are tightly bound in a very compact system and have a rich spectroscopy of excited states. 

We expect three major contributions to the B c decay width: b — ► cW + with the c as 
a spectator, leading to final states like J/tpn or J/ipiv; c —>■ sW + , with the b as specta- 
tor, leading to final states like B s ir or B s £v\ and cb — > W + annihilation, leading to final 
states like D K, r u t or multiple pions. Since these processes lead to different final states, 
their amplitudes do not interfere. When phase space and other effects are included, the 
predicted lifetime is in the range 0.4-1.4 ps IHH. Because of the wide range of predic- 
tions, a B c lifetime measurement is a test of the different assumptions made in the various 
calculations. Several authors have also calculated the B c partial widths to semileptonic final 
states 0-00. 

The production of B c mesons has been calculated in perturbative QCD. At transverse 
momenta large compared to the B c mass the dominant process is that in which a b is 
produced by gluon fusion in the hard collision and fragmentation provides the c. At lower 
Pt a full oQ calculation shows that the dominant process is one in which both the b and 
c quarks are produced in the hard scattering. These calculations [14]-|i~8|| provide inclusive 
production cross sections along with distributions in px and other kinematic variables. 

Experimental searches at LEP have obtained limits on B c production []I9"|-|2~1|| and a 



few candidate events [20,21]. A prior CDF search placed a limit on B c production in the 



B+ — > J/ipir + decay mode |32 

We report here the observation of B c mesons produced in 1.8 TeV pp collisions at the 
Fermilab Tevatron collider using a 110 pb -1 data sample collected with the CDF detector. 
A more detailed description of this work can be found in Ref. ||23|| . We searched for the 
decay channels £>+ — ► J/tp^v and £>+ — ► J/ipe+v followed by J/ip — > fi + fi~. A Monte 
Carlo calculation of B c production and decay to J/ipiv showed that, for an assumed B c 
mass of 6.27 GeV/c 2 , 93% of the J/ip£ final state particles would have J/ip£ masses with 
4.0 < M(J/ip£) < 6.0 GeV/c 2 . We refer to this as the signal region, but we accepted 
candidates with M(J/ip£) between 3.35 and 11 GeV/c 2 . 

We have described the CDF detector in detail elsewhere P^ffiSj] . The tracking system 
for CDF gives a transverse momentum resolution Spx/pr — [(0.0009 x p T ) 2 + (0.0066) 2 ] 1 / 2 , 
where px is in units of GeV/c. The average track impact parameter resolution relative to the 
beam axis is (13 + (40/pr)) /xm in the plane transverse to the beam |26|. An online di-muon 
trigger and subsequent offline selection yielded a sample of about 196,000 J/ip — > n + fi~ 
mesons. 

The decays we sought, B+ — ► J/ip £ + u, have a very simple topology: a decay point for 
J/ip —>■ fi + H~~ displaced from the primary interaction point and a third track emerging from 
the same decay point. This J/tj) + track sample included B+ — > J/ip e + u, B+ — > J/ijj fi + u, 
B + — > J/ipK + , and background from various sources. We subjected the three tracks to a 
fit that constrained the two muons to the J/ip mass and that constrained all three tracks to 
orginate from a common point. A measure of the time between production and decay of a 
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B c candidate is the quantity ct*, denned as 



^_ M(J/j;£).L xy (J/i;£) 

\pt{JM\ [) 

where L xy is the distance between the beam centroid and the decay point of the B c can- 
didate in the plane perpendicular to the beam direction and projected along the direction 
of the J/ipl combination in that plane, M(J/if)£) is the mass of the tri-lepton system, 
and pt(J/iP£) is its momentum transverse to the beam. Our average uncertainty in the 
measurement of ct* is 25 /mi. We required ct* > 60 jum. 

B + — > J/ipK + candidates were identified by a peak in the fi + fi~K + mass distribution 
centered at M(B + ) = 5.279 GeV/c 2 with an r.m.s. width of 14 MeV/c 2 . (See Fig. 2 of 



Ref. ||23|| .) The peak contained 290 ± 19 events after correction for background. Events 
within 50 MeV of M(B + ) were eliminated as B+ — > J/ip£ + v candidates. 

Electrons were identifed by the association of a charged-particle track with p? > 2 GeV/ c 
and an electromagnetic shower in the calorimeter. Additional information for identifying 
electrons was obtained from specific ionization in the tracking chambers and from the shower 
profile in proportional chambers embedded in the electromagnetic calorimeter. Muons from 
J/i/j decay were identifed by matching a charged-particle track with pt > 2 GeV/c to a track 
segment in muon drift chambers outside the central calorimeter (5 to 9 interaction lengths 
thick depending on angle). The third muon was required to have a transverse momentum 
exceeding 3GeV/c and to pass through an additional three interaction lengths of steel to 
produce a track segment in another set of drift chambers. We found 23 B+ — > J/ip e + v 
candidates of which 19 were in the signal region, and we found 14 B+ —>■ J/ip fi + v candidates 
of which 12 were in the signal region. 

Significant contributions to backgrounds in the sample of B c candidates come from 
misidentification of hadron tracks as leptons (i.e. false leptons) and from random combina- 
tions of real leptons with J/ip mesons. There are three significant sources of false leptons: 
hadrons that reach the muon detectors without being absorbed; hadrons that decay in flight 
into a muon in advance of entering the muon detectors; and hadrons that are falsely iden- 
tified as electrons. In one type of random combination external or internal conversions, i.e. 
electrons from photons that convert to e + e~ pairs in the material around the beam line or 
from Dalitz decay of tt°. When the other member of the pair remains undetected, electrons 
from these sources contribute a "conversion background" to the candidate sample. The 
other type of random combination involves a B that has decayed into a J/ip and an asso- 
ciated B that has decayed semileptonically (or through semileptonic decays of its daughter 
hadrons) into a muon or an electron. The displaced J/ip and the lepton can accidentally 



appear to originate from a common point. A number of other backgrounds [23[] were found 



to be negligible. From a combination of data and Monte Carlo calculations, we determined 
the J/ip + track mass distribution for each of the sources of background. As a check of 
our background calculations, we verified that we are able to predict the number of events 
and mass distribution in an independent, background-rich sample of same-charge, low- mass 
lepton pairs. (See Fig. 27 in Ref. [p3f| .) 

We verified the topology for candidate events by applying all selection criteria except the 
requirement that the third track intersect the J / ip vertex. The impact parameter distribution 
between the third track and the J/ip vertex has a prominent peak at zero, demonstrating 
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that, for most candidate events, the three tracks arise from a common vertex. (See Fig. 28 
in Ref. ||.) 

Table | summarizes the results of the background calculation and of a simultaneous fit for 
the muon and electron channels to the mass spectrum over the region between 3.35 and 11 
GeV/c 2 [|23J. Figure |I] shows the mass spectra for the combined J/ipe and J /if) \x candidate 
samples, the combined backgrounds and the fitted contribution from B£ — ► J/ip£ + v decay. 
The fitted number of B c events is 20.4l 5 ; 5 . 

To test the significance of this result, we generated a number of Monte Carlo trials with 
the statistical properties of the backgrounds, but with no contribution from I? c mesons. These 
were subjected to the same fitting procedure to determine contributions consistent with the 
signal distribution arising from background fluctuations. The probability of obtaining a 
yield of 20.4 events or more is 0.63 x 10~ 6 , equivalent to a 4.8 standard-deviation effect. 

To check the stability of the B c signal, we varied the value assumed for the B c mass. 
We generated signal templates, i.e. Monte Carlo samples of B+ — > J/if)£ + v, with various 
values of M(B C ) from 5.52 to 7.52 GeV/c 2 . The signal template for each value of M(B C ) and 
the background mass distributions were used to fit the mass spectrum for the data. This 
study established that the magnitude of the B c signal is stable over the range of theoretical 
predictions for M(B C ), and the dependence of the log-likelihood function on mass yielded 
M(B C ) = 6.40 ± 0.39 (stat.) ± 0.13 (syst.) GeV/c 2 . 

We obtained the mean proper decay length cr and hence the lifetime r of the B c meson 
from the distribution of ct*. We used only events with 4.0 < M(J/if>£) < 6.0 GeV/c 2 , 
and we changed the threshold requirement on ct* from ct* > 60 fim to ct* > —100 //m 
for this lifetime measurement. This yielded a sample of 71 events, 42 J/ipe and 29 J/if> //. 
We determined a functional form for the shapes in ct* for each of the backgrounds. To 
these, we added a resolution-smeared exponential decay distribution for a B c contribution, 
parametrized by its mean decay length cr. Because of the missing neutrino, the proper 
decay length ct for each event differs from ct* of Eq. [I]. We convoluted the exponential in ct 
with the distribution of ct* jet derived from Monte Carlo studies. Finally, we incorporated 
the data from each of the candidate events in an unbinned likelihood fit to determine the 
best-fit value of cr. The data and the signal and background distributions are shown in 
Fig. ||, and the result is: 

cr = 137±|| (stat.) ± 9 (syst) fim (2) 



t = 0.46±g;i| (stat.) ± 0.03 (syst.) ps (3) 

From the 20.4 B c events and the 290 B + — > J/ipK + events, we calculated the B c 
production cross section times the B+ — > J /if) £ + v branching fraction a-BR(B+ — > J '/if) £ + u) 
relative to that for the topologically similar decay B + — > J/if> K + . Systematic uncertainties 
arising from the luminosity, from the J/if> trigger efficiency, and from the track- finding 
efficiencies cancel in the ratio. Our Monte Carlo calculations yielded the values for the 
efficiencies that do not cancel. The detection efficiency for B+ — > J /if) £ + v depends on cr 
because of the requirement that ct* > 60 fim, and we quote a separate systematic uncertainty 
because of the lifetime uncertainty. We assumed that the branching fraction is the same for 
B+ — > J /if) e + v and B+ — > J/ip fi + v. We multiply the 20.4 events by a factor 0.85 ± 0.15 



to correct for other B c decay channels such as B c — > if>' £u |23]. We find 
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= T{B) '■ Tr^Zj'/Ik) = °- 132 -ow (stat.) ± 0.031 (syst.) ™ (lifetime); 

(4a) 

for and B + with transverse momenta pt > 6.0 GeV/c and rapidities \y\ < 1.0. This result 
is consistent with limits from previous searches |T9|-pT|. Figure ^| compares phenomenolog- 
ical predictions with our measurements of cr and TZ(J/ip iv). Within experimental and 
theoretical uncertainties, they are consistent. 

In conclusion, we report the observation of B c mesons through their semileptonic de- 
cay modes, B c — > J/tp£X where I is either an electron or a muon. We measured the B c 
mass and the product of its production cross section times semileptonic branching fraction 
which confirm phenomenological expectations. We measured a B c lifetime consistent with 
calculations in which the decay width is dominated by the decay of the charm quark. 
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Council of Canada; the National Science Council of the Republic of China; the A. P. Sloan 
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TABLES 



TABLE I. B c Signal and Background Summary 



3.35 < M(J/^j£) < 11.0 GeV/c 2 




J/ijie Events 


J ftp n Events 


False Electrons 


4.2 ±0.4 




Undetected Conversions 


2.1 ± 1.7 




False Muons 




11.4 ±2.4 


BB bkg. 


2.3 ± 0.9 


1 44 ± 25 


Total Background (predicted) 


8.6 ± 2.0 


12.8 ±2.4 


(from fit) 


9.2 ±2.0 


10.6 ±2.3 


Predicted N{B C -► J/ip e u)/N{B c -► Jytyfi/) 




0.58 ±0.04 


e and // Signal (derived from fit) 


12.0^;" 




Total Signal (fitted parameter) 




20.4l|i 


Signal + Background a 


21.2 ±4.3 


19.0 ±3.5 


Candidates 


23 


14 


P(Null) b 




0.63 x 10" 6 



a The total number of fitted events was not constrained to be equal to the number of candidates. 
b Probability that background alone can fluctuate to produce an apparent signal of 20.4 events or 
more, based on simulation of statistical fluctuations. 
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FIG. 1. Histogram of the J/ip£ mass that compares the signal and background contributions 
determined in the likelihood fit to the combined data for J/ipe and J/tp^i. Note that the mass 
bins vary in width. The total B c contribution is 20.4^ g events. The inset shows the behavior of 
the log-likelihood function — 21n(L) vs. the number of B c mesons. 
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0. 400. 800. 1200. 

ct Cum) 

FIG. 2. The distribution in ct* for the combined J/ipfi and J/ipe data along with the fitted 
curve and contributions to it from signal and background. The inset shows the logdikelihood 
function vs. ct for the B c . 
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FIG. 3. The point with 1-standard-deviation contour shows our measured value of the a ■ BR 
ratio plotted at the value we measure for the B c lifetime. The shaded region represents theoret- 
ical predictions and their uncertainty corridors for two different values of the semileptonic width 
T s i = T(B C —* J/if)£v) based on Refs. 0] and 0. The other numbers assumed in the theoretical 



predictions are V cb = 0.041 ± 0.005 ||, a(B+)/a(b) = 1.3 x 10~ 3 @, = 0.378 ±0.022 fH, 

-» J/ipK+) = (1.01 ± 0.14) x 10~ 3 
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